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The parton and hadron cascade model PACIAE based on PYTHIA is utilized to systematically 
investigate strange particle production in pp collisions at the RHIC and LHC energies. Globally 
speaking, the PACIAE results of the strange particle rapidity density at mid-rapidity and the 
transverse momentum distribution are better than PYTHIA (default) in comparing with STAR 
and ALICE experimental data. This may represent the importance of the parton and hadron 
rescatterings, as well as the reduction mechanism of strange quark suppression, added in the 
PACIAE model. The K/-n ratios as a function of reaction energy in pp collisions from SPS to LHC 
energies are also analyzed in this paper. 
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I. INTRODUCTION 

Strange particle production is a powerful probe into 
the hadronic interaction and the hadronization process in 
pp and heavy ion collisions at relativistic energies. The 
strangeness enhancement in relativistic nucleus-nucleus 
collisions relative to pp collisions at the same energy has 
been proposed as a signature of the Quark-Gluon Plasma 
(QGP) formation in the relativistic heavy ion collisions 
[l| . This is based on the principle that the threshold en- 
ergy of strange particle production in Quark-Gluon Mat- 
ter (QGM) is higher than that in hadronic matter. Un- 
fortunately, the QGM has not yet been confirmed to be 
unique explanation of strangeness enhancement. 

The Large Hadron CoUider (LHC) at CERN opens a 
new era in the investigation of nucleus- nucleus collisions 
at relativistic energy. It has established a new platform 
to study the properties of QCD matter (QGM) 0, [1]. 
Strangeness production is one of the most important re- 
search topics in the relativistic nuclear collisions at the 
LHC. Recently, the ALICE collaboration reported their 
data of strange particle production in pp collisions at 
= 0.9 TeV [4151. This introduces a new energy regime 
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in the research of strangeness production and the possi- 
bility to compare with the previous measurements in pp 
collisions at Vs=200 GeV 0. 

In the LUND string fragmentation scheme Q , the sup- 
pression of s quark pair production compared with u (d) 
pair production (the parameter parj{2) in PYTHIA or 
A denoted later) was assumed to be fixed. However, 
later experiments f§\ have shown that this suppression 
decreases with increasing reaction energy. In Ref. Q, 
the reduction mechanism of the strange quark suppres- 
sion has been introduced in the LUCIAE model. By this 
mechanism, they described the strangeness enhancement 
in pp, p-l-A, and A-l-A collisions at SPS energies success- 
fully LUCIAE is a hadron cascade model based 
on FRITIOF [ll| with the firecracker model and hadronic 
rescattering added [T2|. 

In this paper, the reduction mechanism of the strange 
quark suppression was introduced in the parton and 
hadron cascade model PACIAE Ts'] and denoted as the 
modified PACIAE model to distinguish it from the de- 
fault one. The modified PACIAE model was then used 
to systematically investigate the strange particle produc- 
tion in pp collisions at the RHIC and LHC energies. 

This paper is organized as follows. In Sec. II, we give 
a brief review the PACIAE model Ts'] and the reduc- 
tion mechanism of strangeness quark suppression Q . In 
Sec. HI, we use the modified parton and hadron cascade 
model PACIAE based on PYTHIA to analyze sys- 
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tematically the strangeness production in pp collisions 
at RHIC and LHC energies [i^. The K/tt ratio as a 
function of reaction energy in pp collisions at relativis- 
tic energies from SPS to LHC was also investigated in 
Sec. III. Section IV gives a summary and conclusion. 



is first created quantum mechanically at one point and 
then tunnel out to the classically allowed region. This 
tunnelling probability is calculated 

exp( )exp( ), (1) 



II. MODELS 

PACIAE Ts^ is a parton and hadron cascade model 
based on PYTHIA [1^] which is a model for high energy 
hadron-hadron (hh) collisions at the hadronic level. In 
the PYTHIA model, a pp collision is decomposed into 
parton-parton collisions. The hard parton-parton inter- 
action is described by the lowest leading order pertur- 
bative QCD (LO-pQCD). The soft parton-parton colli- 
sion, a non-perturbative phenomenon, is considered em- 
pirically. The initial- and final-state QCD radiations, 
as well as the multiparton interactions, are considered. 
Therefore the consequence of a pp collision is a parton 
multijet configuration composed of quarks (anti-quarks), 
di-quarks (anti-diquarks) and gluons, along with a few 
hadronic remnants. After that, the string construction 
and fragmentation (hadronization) are preformed to ob- 
tain a hadronic final state for a pp collision. 

For the pp collision, the PACIAE model is different 
from PYTHIA in the addition of the parton initiation 
stage, the parton rescattering before hadronization, and 
the hadron rescattering after hadronization. Thus, the 
PACIAE model consists of the parton initiation, parton 
evolution (rescattering), hadronization, and hadron evo- 
lution (rescattering) four stages. 

In order to create the parton initiation stage for a pp 
collision, the string fragmentation is switched off tem- 
porarily in PACIAE. The di-quarks (anti-diquarks) are 
broken up. Then the partonic initial state is obtained in 
the PACIAE model, instead of the hadronic final state in 
PYTHIA. This is just QGM formed in the parton initia- 
tion stage of a pp collision. 

The rescattering among partons in QGM is considered 
by the 2 — >■ 2 LO pQCD parton-parton cross sections [l5| . 
By integrating this differential cross section properly, the 
total cross section is obtained. With the differential and 
total cross sections, parton rescattering is performed by 
the Monte Carlo method. 

In the hadronization stage, the QGM formed after par- 
ton rescattering is hadronized by the LUND string frag- 
mentation regime [7^14| or the Monte Carlo coalescence 
model proposed in flsfT^The LUND string fragmentation 
is used in this paper. 

The hadronic rescattering is modeled with the 
usual two body elastic and/or inelastic collisions 
[T^ . until the hh collision pairs are exhausted 
(hadronic freeze-out). The rescatterings among 
TT, K, p, n, p(u}), A, A, E, 5, J/^P and their antipar- 
ticles are considered for the moment. 

In the PYTHIA model [11], it is assumed that the qq 
pair with quark mass m and transverse momentum pT 



where the string tension is assumed to be a constant of 
K «1 GeV/fm« 0.2 GeV^ 0,111. This probability im- 
plies a suppression of strange (heavy) quark production: 
u : d : s : c 1 : 1 : 0.3 : 10^^^. Therefore the charm 
and heavier quarks are not expected to be produced in 
the soft string fragmentation process but only in the hard 
process or as a part of the initial- and final-state QCD 
radiations. 

A reduction mechanism of strange quark suppression 
was introduced in [§] by assuming that the effective string 
tension increased with increasing reaction energy. Hence 
the strange (heavy) quark production increased with in- 
creasing reaction energy. It was further assumed in Q 
that the effective string tension variation with reaction 
energy could be considered by the effective string tension 
as a function of the number and the hardening of gluons 
in a single string as follows: 



^^0(1 -6" 
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where kq is the string tension of a pure qq string, as- 
sumed to be ~ 1 GeV/fm. Here it should be mentioned 
that the above Eq. ([3]) represents the deviation scale of 
the multigiuon string from the pure string. The gluons in 
a multigiuon string are ordered from 2 to rt — 1 , because 
of the quark and antiquark on both ends of string with 
index 1 and n, respectively, k^^ (kTmax) is the trans- 
verse momentum of gluon j with kxj > sq (gluon largest 
transverse momentum). The parameters a=3.5 GeV and 
y^=0.8 GeV were determined by fitting the hh collision 
data 0. 

The strange quark suppression factor A (i.e. parj(2) 
in PYTHIA) and the width a {par j (21) in PYTHIA) 
of the Gaussian p^ and Py transverse momentum dis- 
tributions of the primary hadrons, in a string with ef- 
fective string tension k^^'^, are denoted by Ai and cri, 
respectively. These two quantities in a string with effec- 
tive string tension k'^'^^ , A2 and (J2, can be calculated by 
Eq. (P 



A, = a: 
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In the PYTHIA model there are parameters parj(l) 
and parj(3) related to strangeness production, be- 
sides parj(2) and parj(21). parj(l) is the suppression 
of diquark-antidiquark pair production relative to the 
quark- ant iquark pair production. parj(3) refers to the 
extra suppression of strange diquark production relative 
to the normal suppression of strange quark production. 
For double strange particle (strange baryon) production, 
parj(l) and parj(3) have to be considered as well. It is 
not hard to prove that the Eq. ^ is also valid for parj{l) 
and parj(3). 



III. CALCULATIONS AND RESULTS 

The reduction mechanism of the strange quark sup- 
pression has been included in the PACIAE model. One 
might first tune the parameters parj{l), parj(2), and 
parj(3) to fit the strangeness production data in a given 
nuclear collision system at a given energy. The result- 
ing parj{l), parj(2), and parj{3) and the effective string 
tension can be used to predict the strangeness production 
in the same reaction system at different energies, even in 
different reaction systems. 

As we aimed at the physics behind the strange parti- 
cle production and not at the reproduction of experimen- 
tal data, we kept all model parameters at their default 
values. However, the higher order and non-perturbative 
correction factor in parton-parton interaction, K, was 
assumed to be 3. 

We then first globally tuned the parameters parj{l), 
parj{2), and parj{3) in default PACIAE simulations to 
fit the strangeness production data from NSD (non sin- 
gle diffractive) pp collisions at ^/s—200 GeV [6]. The 
results are shown in Tab. U along with the results from 
default PYTHIA simulations. From this table, we find 
that the STAR data are better described by default PA- 
CIAE model rather than by the default PYTHIA model. 

The strange particle transverse momentum spectra in 
the NSD pp collisions at y^=200 GeV calculated by 
the default PACIAE and default PYTHIA are shown in 
Fig. [TJ The corresponding STAR data [1] are also given. 
The panels (a), (b), and (c) are for K°, A, and S~, re- 
spectively. These results indicate that, globally speaking, 
the strange particle pr spectrum is better described by 
the default PACIAE rather than the default PYTHIA 
model. Here, the role of parton and hadron rescatterings 
is also not negligible. 

The fitted parameters of par 0.18, parj{2)—0.4:, 
and parj{3)=0.5 as well as the calculated k'^^-^ =1.387 
were then used to calculate the strange particle produc- 
tion in INEL (inelastic) pp collisions at y/s=900 GeV 
by the modified PACIAE model. The resulting strange 
particle rapidity densities at mid-rapidity were compared 
with the ALICE data [1] as well as the results of default 
PYTHIA in Tab. HH One sees in this table that the AL- 
ICE data are better reproduced by modified PACIAE, 
rather than the default PYTHIA model. The effects of 



the reduction mechanism of strange quark suppression 
and the parton and hadron rescatterings may be all im- 
portant which has to be separately studied. 

We present the strange particle transverse momentum 
distributions in INEL pp collisions at ^/s=900 GeV cal- 
culated by the modified PACIAE and default PYTHIA 
models in Fig. [2] The corresponding ALICE data Q 
are also given. The panels (a), (b), (c), and (d) are for 

K^, (j), A, and S^+S^, respectively. One sees in Fig. [5] 
that the modified PACIAE results fit well with the AL- 
ICE data for the light strange particle if". However, 
for the meson and the heavy strange baryons (A and 

S~-|-S''') the modified PACIAE model overestimates the 
data in the pt >1 GeV/c region. This may be because 
we only used the default PACIAE to tune the parameters 
to the data of strange particle rapidity densities at mid- 
rapidity in -^3=200 GeV pp collisions [6,]. In this sense, 
the transverse momentum distribution is calculated pa- 
rameter free. Therefore, even if the calculated strange 
particle rapidity densities at mid-rapidity are close to 
the data (cf. Tab. H] and Tab. |n|, the calculated trans- 
verse momentum distribution may not be closely similar 
to the data. Thus, the calculated transverse momen- 
tum distribution is softer than the data in the pt <1 
GeV/c region, but harder in the pr >1 GeV/c region. 
It might also be possible that the reduction mechanism 
of strange quark suppression is not so suitable for the 
heavy strange baryons. Therefore, the role played by the 
reduction mechanism of strange quark suppression and 
the parton and hadron rescatterings are worthy to be 
separately studied in the next investigations. 

Similarly, the transverse momentum distribution of 
ifj. A, and in INEL pp collisions at •\/s=2.36, 7, and 
14 TeV were predicted by the modified PACIAE model 
and given in Fig. |31 These results are prepared to com- 
pare with the future experimental measurements at the 
LHC. 

In [5', Tg] , the experimental kaon to pion ratios in pp 
and/or pp collisions at different reaction energies from 
SPS to LHC were plotted in a K/n verses figure, 
although their experimental details in the measured ob- 
servable {K+/n+, or {K+ + K~)/{tt+ + TT-), or K^/n") 
and the rapidity acceptance were different from each 
other. These experimental ratios are copied in Fig. 2] In 
this figure, NA49, STAR ALICE, E735, and UA5 data 
were taken from (H] , [l|| , [1, [11] , [li| , and [131 , respec- 
tively. In @, [ill, it was concluded that the "ratio shows a 
slight increase from y/^=200 GeV {K/n = 0.103 ±0.008) 
to ^=900 GeV {K/n = 0.123 ± 0.004 ± 0.010), yet 
consistent within the error bars. The results at 7 TeV 
will show whether the K/n ratio keeps rising slowly as a 
function of ^/s or saturates." To respond to this chal- 
lenge the (7^+ + K~)/{tt^ -I- 7r~) ratio was extracted 
from the above PACIAE model simulations for pp col- 
lisions at Vs=0.2, 0.9, 2.36, 7, and 14 TeV and given 
in Fig. |4] by red triangles. In addition, the default PA- 
CIAE and PYTHIA models were used to calculate the 
NA49 {K+ +K^)/{'K+ +TT~) ratio in pp collisions at 158 
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GeV/c beam momentum. The model parameters em- 
ployed in these calculations were the same as the ones 
used in the above simulations, except K=l was assumed 
because of the low energy. 

We find in Fig.Hthat the PACIAE results at ^^=17.2, 
200 and 900 GeV agree with the NA49, STAR, and AL- 
ICE data, respectively, within error. Together with the 
PACIAE resuhs at v^=2.36, 7, and 14 TeV, it seems 
to indicate that the K/ir ratio as a function of ^/s first 
increases slightly from -^5=0. 2 to 0.9 TeV and then satu- 
rates. Of course, this rough trend of the /n^ ratio as 
a function of ^/s in pp reactions has to be further studied 
in detail both experimentally and theoretically. 

IV. CONCLUSIONS 

In summary, we used the parton and hadron cascade 
model PACIAE with the reduction mechanism of the 
strange quark suppression and default PYTHIA models 
to systematically investigate the production of strange 
particles in NSD and INEL pp colhsions at y^=200 GeV 
and y/s=0.9, 2.36, 7, and 14 TeV, respectively. The cal- 
culated strangeness transverse momentum distributions 
in NSD pp collisions at y/s=20Q GeV compare well with 
the STAR data ^ and PACIAE is globally better than 
the default PYTHIA model. We also compared the cal- 
culated strange particle rapidity densities at mid-rapidity 
and transverse momentum distributions in INEL pp col- 
lisions at TeV with the corresponding ALICE 
data 0]. The agreement between theory and experiment 
is good and the PACIAE is globally better than the de- 
fault PYTHIA model. This may indicate the significant 
effect of the parton and hadron rescatterings, as well as 
the reduction mechanism of strange quark suppression. 
These effects are worth investigating separately. 

The K/n ratios calculated by PACIAE in NSD pp colli- 



sions at y/s=Q.2 TeV and in INEL pp collisions at y/s=0.9 
TeV agree with the corresponding STAR [3| and AL- 
ICE data respectively, within error. The PACIAE 
model result in pp collisions at 158 GeV/c beam momen- 
tum is also comparable with NA49 data fl7|. The rough 
trend of the K/tt ratio as function of ^/s shown in Fig. 0] 
seems to be first increasing moderately, then slightly, and 
finally saturating. This trend has to be studied in detail 
both experimentally and theoretically. 

Just before the submission of this paper, we read the 
paper of [2lj and knew that the CMS measured "pt dis- 
tributions are found to differ substantially from PYTHIA 
results and the production rates exceed the predictions 
by up to a factor of three." Note, here PYTHIA refers 
to a variety of PYTHIA tunes quoted in [2lj. We use 
the modified PACIAE model to calculate the correspond- 
ing values. The calculated dN/dy\y~o= 0.238, 0.102, 
and 0.0113 for K^, A, and S~ compare well with the 
CMS data of 0.205 ±0.001 ±0.015, 0.108 ±0.001 ±0.012, 
0.011 ±0.001 ±0.001 [Hi in NSD pp collisions at ^/i=0.9 
TeV, respectively. Similarly, the PACIAE results of 
0.403, 0.166, and 0.0183 compare well with the data 
of 0.346 ± 0.001 ± 0.025, 0.189 ± 0.001 ± 0.022, and 
0.021±0.001±0.003, respectively, in NSD pp collisions at 
y/s=7 TeV. The calculated rapidity and transverse mo- 
mentum distributions are compared with the CMS data 
(21} in Fig. [SI The agreement is better than the PYTHIA 
tunes quoted in j2 Ij . 
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TABLE I: Strange particles rapidity densities at mid-rapidity {\y\ <0.5) in NSD pp collisions at y^=200 GeV. The STAR data 
were taken from :6]. 





STAR 


PACIAE 


PYTHIA 




0.134 ±0.011 


0.127 


0.107 


K+ 


0.140 ±0.01 


0.135 


0.112 


R- 


0.137 ±0.009 


0.123 


0.104 


A 


0.0385 ± 0.0036 


0.0360 


0.0163 


A 


0.0351 ± 0.0033 


0.0350 


0.0163 




0.0026 ± 0.00092 


0.00373 


0.00106 




0.0029 ± 0.001 


0.00364 


0.00099 


- +1^ 


0.00034 ± 0.00019 


0.00026 


0.00005 



TABLE II: Strange particles rapidity densities at mid-rapidity in INEL pp collisions at -^8=900 GeV. The ALICE data were 
taken from 4]. 



\y\ pT{GeV/c) ALICE PACIAE PYTHIA 

K^s < 0.75 [0.2-3.0] 0.184 ± 0.0063 (U5 0.172 

$ < 0.60 [0.7-3.0] 0.021 ±0.005 0.019 0.0134 

A < 0.75 [0.6-3.5] 0.048 ± 0.00412 0.043 0.0227 

A < 0.75 [0.6-3.5] 0.047 ± 0.00539 0.043 0.0230 

H-±H"^ < 0.80 [0.6-3.0] 0.0101 ± 0.00219 0.0086 0.00354 
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FIG. 1: (Color online) Transverse momentum distribution of 
strange particles in NSD pp collisions at -^3=200 GeV. Panels 
(a), (b) and (c) are for , A, and H^, respectively. The 
STAR data were taken from [a|. 
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FIG. 2: (Color online) Transverse momentum distributions of strange particles in INEL pp collisions at y^=900 GeV. Panels 
(a), (b), (c), and (d) are for K^, cf>, A, and S~+H , respectively. The ALICE data were taken from 
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FIG. 3: (Color online) Transverse momentum distributions of 
strange particles in INEL pp collisions at -^8=2. 36, 7, and 14 
TeV calculated by PACIAE. Panels (a), (b), and (c) are for 
K° {\y\ <0.75), A {\y\ <0.75), and H" <0.80), respec- 
tively. 
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FIG. 4: (Color online) K/n ratio as a function of ^/s in pp 
collisions at relativistic energies. 
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FIG. 5: (Color online) Left panel: the strange particle rapidity distributions for (a), A (b), and H (c) in NSD pp collisions 
at -^3=0. 9 and 7 TeV. Right panel: the transverse momentum distributions for (d), A (e), and H~ (f). The CMS data were 
taken from f2l|. 



